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An analysis of the dilepton measurements in the reactions Pb(158 A-GeV) + Au, Pb 
points to a unique thermal source contributing to the invariant mass and transverse mo- 
. mentum spectra. Effects of the flow pattern are discussed. 

in 



> 
\o 

(N 
oo . 

O . Introduction: Dileptons are penetrating probes which carry nearly undisturbed in- 

^ I formation about early, hot and dense matter stages in relativistic heavy-ion collisions. 
p ■ Some effort, however, is needed for disentangling the various sources contributing to the 
Qh! observed yields and for identifying the messengers from primordial states of strongly in- 
^ • teracting matter. 

- ^ , The dielectron spectra for the reaction Pb(158 A-GeV) + Au measured by the CERES 

^ • collaboration cannot be described by a superposition of e~^e~ decay channels of final 
hadrons, i.e. the hadronic cocktail. A significant additional source of dielectrons must 
be there. Since the data |jl| cover mainly the invariant mass range M < 1.5 GeV the 
downward extrapolation of the Drell-Yan process is not an appropriate explanation. Also 
correlated semileptonic decays of open charm mesons have been excluded . As a widely 
accepted explanation, a thermal source is found to account for the data (cf. |] and 
further references therein, in particular with respect to in-medium effects and chiral sym- 
metry restoration). 

Very similar, the NA50 collaboration has found, for the reaction Pb(158 A- GeV) + 
Pb, that the superposition of Drell-Yan dimuons and open charm decays does not explain 
the data in the invariant mass range 1.5 GeV < M < 2.5 GeV [^]. Final state interactions 
PI or abnormal charm enhancement |^ have been proposed as possible explanations. 
Here we try to explain the NA50 measurements by a more apparent idea ^, namely 
a thermal source. We present a schematic model which at the same time describes the 



CERES and NA50 data. 

The model: Since we include the corresponding detector acceptances a good starting 

point for Monte Carlo simulations is the differential dilepton spectrum 
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where Q = P1+P2 is the pair four-momentum, pi 2 are the individual lepton four-momenta 

composed of transverse momenta p± 1^2 and rapidities yi^2 and azimuthal angles 0i,2- Here 

we extensively employ the quark - hadron duality 0, ^] and base the rate R on the 

lowest-order quark - antiquark (qq) annihilation rate (cf. [|TU], ^) 

dR ( u ■ Q'] , , 

where u{x) is the four- velocity of the medium depending on the space-time as also the 
temperature T{x) does. Note that, due to Lorentz invariance, u necessarily enters this 
expression. The above rate is in Boltzmann approximation, and a term related to the 
chemical potential is suppressed. As shown in the qq rate deviates from the hadronic 
one at M < 300 MeV, but in this range the Dalitz decays dominate anyhow; in addition, 
in this range the thermal yield is strongly suppressed by the CERES acceptance. In the 
kinematical regions we consider below, the lepton masses can be neglected. 

Performing the space-time and momentum integrations in eqs. H) one gets 
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, 47r f drr"^ for 3L', 
V{t) = { ' ' (5) 

t J drr for 2D, 

where V{t) acts on E, and 3D means spherical symmetric expansion, while 2D denotes 
the case of longitudinal boost-invariant and cylinder-symmetrical transverse expansion; 



the quantity ^ is defined as ^ = T ^ sinh pJ M]_ coshy^ — M^, and p(r, t) is the radial or 



transverse expansion rapidity; Kq and Iq are Bessel functions [|T0[ . The components of the 
lepton pair four-momentum Q = {M± cosh Y, M± sinh Y, Q±) are related to the individual 
lepton momenta via 

Ml = pii +PI2 + 2p±iP±2Cosh(?/i - 7/2), (6) 

Q± = P±l+P±2, (7) 

= Ml-Ql, (8) 
tanhF = P^i sinh|/i+p^2 sinhy2 _ .gx 
p±i cosh yi + p^2 cosh y2 ' 



It turns out that the shape of the invariant mass spectrum dN / [dM dY\\Y\^Q^^dtdV{t)), 
which is determined only by the emissivity function E, does not depend on the flow 
rapidity p in the 2D case |T^, and in the 3D case for T = 120 ■ ■ ■ 220 MeV and p < 0.6 
there is also no effect of the flow. The analysis of transverse momentum spectra of various 



hadrons species point to an average transverse expansion velocity v± ^ 0.43 at kinetic 
freeze-out, while a combined analysis of hadron spectra including HBT data yields v± ~ 



0.55 Therefore, p < 0.6 is the relevant range for the considered reactions. 

We note further that the invariant mass spectra dN/ {dM dY\\Y\<o.5 dt dV{t)) for the 
3D and 2D cases differ only marginally. Relying on these findings one can approximate 
the emissivity function E by that of a "static" source at midrapidity, as appropriate only 
for symmetric systems, 

f M^coshFl , , 

thus getting rid of the peculiarities of the flow pattern. In contrast to the invariant 
mass spectra, the transverse momentum or transverse mass spectra are sensitive to the 
flow pattern |Ty, 0, |T3], in general. A value of p = 0.4, for example, causes already 
a sizeable change of the shape of the spectra dN/{dQ±dY\\y^^Q^5dtdV(t)) compared to 
p = 0, in particular in the large-Q^ region. The differences between the 2D and 3D cases 
are not larger than a factor of 2 and, in a restricted Q± interval, can be absorbed in a 
renormalization. The most striking difference of the 2D and 3D scenarios is seen in the 
rapidity spectrum: for 2D it is flat, while in the 3D case it is localized at midrapidity 
(values of p < 0.6 also do not change the latter rapidity distribution). Below we shall 
discuss which space is left to extract from the dilepton data in restricted acceptance 
regions hints for the flow pattern when the other dilepton sources are also taken into 
account. 

Comparison with data: In line with the above arguments we base our rate calcu- 



lations on eqs. (y, |T^) and use the parameterizations ||T6| 

T = (T,-roo)exp|-^|+Too, (11) 

V = A^expj^j. (12) 

with Ti = 210 MeV, = 110 MeV, ti = 5 fm/c, ta = 8 fm/c, N = |±g with A,B as 
mass numbers of the colliding systems and no = 0.17 fm~^. We stop the time evolution 
at Tf = 130 MeV. 

In fig. 1 we show the comparison with the preliminary CERES data applying the 
appropriate acceptance |jl|]. One observes a satisfactory overall agreement of the sum 
of the hadronic cocktail and the thermal contribution with the data. It is the thermal 



contribution which fills the hole of the cocktail around M = 500 MeV in the invariant 
mass distribution in fig. la. In the mass bin M = 0.25 ■ ■ • 0.68 GeV the thermal yield is 
seen (fig. lb) to dominate at small values of the transverse momentum Q±. In this region 
of Q± transverse fiow effects are not important. The large-Q_L spectrum is dominated by 
the cocktail. For higher-mass bins the thermal yield in the region of the first two data 
points is nearly as strong as the cocktail and rapidly falls then at larger values of Q± 
below the cocktail. Therefore, the fiow effects turn out to be of minor importance for the 
present analysis, since within our framework the transverse flow shows up at larger values 
of Q^. 

The question now is whether the same thermal source model accounts also for the 
NA50 data ||^. The available NA50 data are not yet efficiency corrected and the efficiency 
matrix is not accessible. To have a reference we proceed as follows. According to the 
detailed analysis ^ of the most central collisions, the shapes of the M and Q± spectra 
for invariant masses below the J/ip peak are described by the Drell-Yan yield + 3 x the 
yield from correlated semileptonic decays of open charm mesons, both ones generated 



with PYTHIA [0]. We have checked that our PYTHIA simulations coincide with results 
of the NA50 collaboration when applying the acceptance cuts. Our K factors are adjusted 
by a comparison with Drell-Yan data [l^ and identified open charm data (cf. for the 
data compilation). Some confidence in our procedure is gained by correctly reproducing 
the NA50 interpretation of the NA38 data []18[, which are efficiency corrected and 
which we can directly analyze: the NA38 data are described by the Drell-Yan yield -|- 
1.45 X the open charm contribution; the absolute normalization is obtained from a fit of 
the Drell-Yan region beyond the J/ip. 

To get the yield for Pb + Pb collisions from PYTHIA the overlap function T^a = 31 
mb~^ is used. The resulting spectra (within the kinematical cuts for the NA50 acceptance) 
are displayed in figs. 2a and 2b by open squares. One observes that both the continuum 
part of the invariant mass spectrum (without the J/ip contribution) and the transverse 
momentum spectrum for the mass bin M = 1.5 ■ ■ ■ 2.5 GeV are nicely reproduced by the 
sum of Drell-Yan, open charm and thermal contributions. The thermal yield dominates 
again at not too large values of Q± where transverse fiow effects can be neglected. There- 
fore, it seems that from present dilepton measurements the transverse flow can hardly be 
inferred. 

Summary and discussion: In summary we have shown that a very simplifled, 
schematic model for thermal dilepton emission, with parameters adjusted to the CERES 
data, also accounts for the measurements of the NA50 collaboration. A direct comparison 



with the NA50 data is not possible as long as the efficiency matrix is not at disposal. 
Nevertheless, our study points to a common thermal source seen in different phase space 
regions in the dielectron and dimuon channels. This unifying interpretation of different 
measurements has to be contrasted with other proposals of explaining the dimuon excess 
in the invariant mass region 1.5 ■ ■ ■ 2.5 GeV either by final state hadronic interactions or 
an abnormally large open charm production. The latter one, in principle, can be checked 
experimentally 0, thus improving our understanding of dilepton sources. 

Due to the convolution of the local matter emissivity and the space-time history of the 
whole matter and the general dependence on the flow pattern, it is difficult to decide which 
type of matter (deconfined or hadron matter) emits really the dileptons. Our model is not 
aimed at answering this question. Instead, with respect to chiral symmetry restoration, 
we apply the quark - hadron duality as a convenient way to roughly describe the dilepton 
emissivity of matter by a qq rate, being aware that higher-order QCD processes change 
this rate IjlR EU| (to some extent this might be included in a changed normalization A^). 



In further investigations a more reliable rate together with a more detailed space-time 
evolution must be attempted and chemical potentials controlling the baryon and pion 
densities must be included. In this line the model has to be improved before using it for 
analyzing the dilepton yields in the asymmetric reactions S(200 A- GeV) + X measured 
by the NA38, HELIOS-3, and CERES collaborations, where also details of the rapidity 
distribution become important. 
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Figure 1: The preliminary CERES data [||] and the hadronic cocktail [|l| (dashed lines) 
and the thermal yield (full curves). The sum of the cocktail and the thermal yield is 
shown by the gray curves. Left panel (a): the invariant mass spectrum, right panel (b): 
the transverse momentum spectrum for the mass bin 0.25 ■ • ■ 0.68 GeV. 




Figure 2: The reconstructed measurements of the NA50 collaboration (see text for 
details) indicated by open squares (not data!) in comparison with the thermal yield (full 
curves), the Drell-Yan contribution (dashed curves) and the contribution of open charm 
decays (dash-dotted curves). The sum of these contributions is displayed by the gray 
curves. Left panel (a): the continuum invariant mass spectrum, right panel (b): the 
transverse momentum spectrum for the mass bin L5 ■ ■ ■ 2.5 GeV. 



